Possible Localized Modes in the Uniform Quantum Heisenberg Chains of Sr2Cu03. 
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A model of mobile bond-defects is tentatively proposed to analyze the "anomalies" observed on the 
NMR spectrum of the quantum Heisenberg chains of Sr2Cu03. A bond-defect is a local change in 
the exchange coupling. It results in a local alternating magnetization (LAM), which, when the defect 
moves, creates a flipping proces of the local field seen by each nuclear spin. At low temperature, 
when the overlap of the LAM becomes large, the defects form a periodic structure, which extends 
over almost all the chains. In that regime, the density of bond-defects decreases linearly with T. 
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In low-dimensional quantum antiferromagnets, the ef- 
fects of impurities or defects in "gapped" and "un- 
gapped" systems may display drastic differences but also 
common features [0. A good example of the latter case 
is provided by the uniform quantum (s = 1/2) Heisen- 
berg chains (UQHC). Hereafter, we consider the effects 
due to defects in such systems. A "defect" is defined as 
a local change in the magnetic bond coupling. As the 
translational symmetry of the spin system is broken, a 
local alternating magnetization (LAM) develops around 
the defects in the direction of the applied field: this is 
one of the common features mentioned above jl] . Such a 
LAM is probed accurately by NMR [|j Our analysis 
relies on results obtained on the compound Sr2Cu03 [Q. 
In such a "real system" , the elastic coupling between the 
spins and the lattice cannot be ignored. It can result in a 
lattice distortion giving rise to the so-called spin-Peierls 
transition (at temperature Tgp) Dynamical effects 
on LAM structures can also been expected from such 
magneto-elastic couplings. In the dimerized (D) phase of 
a spin-Peierls system (SPS), i.e., for T < Tsp, the LAM 
induced by chain-end effects become mobile in the pure 
one-dimensional case Q. In the "high-field" phase of 
a SPS, a periodic LAM structure occurs, which results 
from the field-induced incommensurate lattice modula- 
tion 1^. In that case, dynamical effects are due to the 
quantum fluctuations of the phason modes describing the 
vibrations of the modulated lattice structure In the 
uniform phase of a SPS, i.e. for T > Tgp, the spin-lattice 
coupling is also of crucial importance. It can apprecia- 
bly change the magnetic susceptibility of the system well 
above Tsp §. 

In Sr2Cu03, no spin-Peierls transition is observed as a 
3-dimensional magnetic ordering takes place at T^r — 5 
K. Above Ttv, this compound provides a remarkable ex- 
ample of UQHC. The exchange coupling is large (J ~ 
2200 K). The logarithmic corrections characterizing the 
quantum (T << J) behavior of UQHC have been ob- 
served M . LAM associated with the edges of finite chains 



have also been observed in Sr2Cu03 Ig]. Due to that 
LAM, the NMR lineshape is changed in a very specific 
way, giving rise to "features" , which have been well iden- 
tified. This is the case of features "A" displayed in Figs. 
1 and 2. Other features, however, are visible on the NMR 
line (features "B" and "C" in Figs. 1 and 2), which have 
not been explained yet. We propose an analysis of these 
additional features in terms of LAM associated with mo- 
bile "bond-defects" . 

The NMR spectra to be analyzed - a few examples 
are displayed in Figs. 1 and 2 - have been obtained on 
an Ar-annealed single crystal. That annealing procedure 
minimizes in that compound the possible source of im- 
purities induced by interstitial excess oxygen. After such 
a treatment, the concentration of the residual spin-1/2 
impurities becomes as low as 1.3 10"'' Some of the 
NMR spectra displayed in Figs. 1 and 2 have been pre- 
viously reported (Figs, la, b and Fig. 2a) in |^ together 
with the experimental conditions. As clearly established 
in 1^, features "A" originate from the chain edges, while 
we show here that features "B" is well explained by the 
presence of mobile bond-defects. Features "C" would re- 
sult from the periodic structure they form at low T. Fi- 
nally, possible origins of these defects are suggested in 
relation with the lattice properties. 

First we consider features "A". LAM associated with 
the edges of finite UQHC have been described by Eg- 
gert and Affleck [|j . The LAM amplitude increases with 
distance from the chain end. At T 7^ however, ther- 
mal fluctuations act as a cutoff, and at long distance, the 
LAM amplitude decreases exponentially with a charac- 
teristic scale given by the correlation length of the sys- 
tem: ^ = J/2T. An analytical expression for the local 
susceptibility has been derived, which agrees well with 
Montc-Carlo calculations. Accordingly, for an hyperfine 
coupling (expressed in GHz and /3 = a, b or c refer- 
ring to the field direction with respect to the crystal axes 
- chain axis b), the local field seen by a nuclear spin at 
site n (counted from a chain end) can be written 
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h,, = -{-irB Ap H {{n + (/.)/ Vesinh(2(n + 0)/^ (1) 

with B ~ 0.020 GHz^^ and /i„ and H expressed in Tesla. 
In (1), (j> allows for a small shift of the LAM along the 
chain, (jj changes the position rimaxisee Fig. 3a) of the 
maximum field value, but for (f) < Umax, it does not af- 
fect this value: h^ax ^ 1-08 10"^ Ap Hy^ {Aa ~ 0.14 
GHz and A^ ~ 0.50 GHz). In the NMR spectrum, the 
field splitting AH (see Fig. 1) where features "A" occur 
is given by AH = 2hmax- In that respect, as explained 
in 1^ , a very good agreement is obtained between theory 
and experiment. The shape of features "A", however, 
depends appreciably on (j) (see dot line = in Fig. lb). 
Using (p as an adjustable parameter, a good agreement 
with the experimental line shape can be achieved for all 
the collected data (see dots in Figs. 1 and 2). With that 
procedure, one obtains that cj) is field independent. It 
shows, however, a T dependence: (j) = (1200 ± 200)/T, 
which, remarkably, agrees with </> ~ ^. As shown in Fig. 
3a (solid line), for that particular value of cf), the LAM 
displays no maximum, but a rather "flat" initial behav- 
ior. Despite the rough description used here, one is led to 
conclude that the initial increase of the LAM amplitude 
is not observed experimentally, in contradiction with the 
theoretical prediction The elastic spin-lattice cou- 
pling may play here an important role. The model of a 
sudden cutoff of the exchange coupling (J = 0) just at 
the chain edges would not applied to such "real" systems. 
Instead, as in the D phase of a SPS [|oj, on a length scale 
~ ^, a distribution of the J-bonds would take place, giv- 
ing rise to the observed "smoothing" of the initial LAM 
amplitude. 

Second, we consider features "B". We assume that 
they result from a local change in the exchange coupling 
(J' 7^ J). As studied by Eggert (for J' < J) LAM 
similar to the case J' = develop around such bond- 
defects. We suppose that the associated LAM can be 
described by an expression similar to (1), but with a 
multiplying factor a (a <1) to account for the expected 
smaller amplitude. As before, a parameter (p' describes 
a possible "smoothing" effect around the defect position. 
Finally, we allow these defects to move along the chains. 
This dynamical behavior is of crucial importance for the 
NMR lineshape as it results in a "flipping" process of the 
hyperfine local field /i^ seen by the nuclear spins. Accord- 
ingly, the NMR signal of each nuclear spin is composed 
of 1 or 2 distinct lines depending on the flipping rate G 
with respect to the local hyperfine field h'^. That "mo- 
tional narrowing" is described as follows. For a nuclear 
spin at site n in the chain, the NMR spectrum is written 

/„ = [Do + K' + K" {u)f +DI+K' [uf] (2) 

with K' [v) = Ylc{K{v)}, K" [v) ^ Iwl{K{v)} and 
V = v'^ — "jH/2'K, where v'^ is the experimental fre- 
quency and Dq is the "natural" NMR width. K{y) de- 
scribes the effect of the dynamical local field as K{v) = 



{'jh[^/2Tr)'^k{i'). Here, k{i') is the Laplace representation 
of k{t) which, as a function of time t, describes the flip- 
ping process due to the motion of the defects. A balhstic 
motion results in a linear exponential, k{t) = exp{—t/G), 
while a diffusive behavior gives k{t) = exp{—^yt/G) [ p2[ . 
Such a single moving defect in a long chain may cre- 
ate additional "features" in the total NMR line, i.e., 
/ = Tin In- That model allows us to reproduce very 
well the observed features "B". An example is given in 
Fig. 3b, where a comparison with the data obtained for 
// a, at i/° = 87 10"^ qh^ and T = 30 K is presented. 
The left side corresponds to a diffusive behavior with the 
parameters: 4>' = 20, a = 0.6, G = 2.8 10"^ GHz and 
Dq = 1.5 10~^ GHz. The ballistic model (right side) gives 
usually a slightly poorer agreement (here with a = 0.45, 
G = 1.5 10"'' GHz). That procedure apphed to the dif- 
ferent experimental conditions yields an interesting result 
concerning 0'. This parameter is independent on H, but 
it depends on T according to cp' = (600 ± 100)/T ~ ^/2. 
As for ip, this smoothing effect can be assumed to result 
from a local induced lattice distortion around the bond- 
defects. 

In presence of several such bond-defects, the above de- 
scription applies as long as the average distance between 
the defects remains large compared to ^, i.e., a model 
of independent defects. At low T, however, interactions 
between defects are expected: a magnetic interaction 
should result from the overlap of the LAM structures and 
an elastic interaction from the lattice distortion associ- 
ated with the defects. As a balance between these two 
interactions, we assume now that the defects form a peri- 
odic structure. If it is the case, the shape of the resulting 
LAM will depend on the number / of spins between two 
neighboring bond-defects. From the dashed line in Fig. 
3a, approximate descriptions of such a LAM are repre- 
sented in Fig. 3c, for I odd and even, respectively. As 
expected, for / even, a node occurs at the middle posi- 
tion. For I odd, the LAM amplitude remains finite at 
all positions. In the former case, the resulting NMR line 
displays a single narrow peak (Fig. 3d, right side). In 
the latter case, however, as the local field seen by the 
nuclear spins never cancels, a double structure develops 
at the centre of the NMR line (left side). For 7J // a , 
in particular (see Fig. 1), such a central double struc- 
ture is clearly observed below 30 K (and well reproduced 
by our model). For H / / c, the NMR line observed at 
T = 30 K (Fig. 2a) is more complex. In addition to 
the overlap of the three components of the quadrupolar 
splitting 1^, it can be analyzed as composed of two con- 
tributions as illustrated in Fig. 2b. For one contribution 
(the dotted line), each component displays a narrow cen- 
tral peak (similar to the case described in Fig. 3d, right, 
for / even) and for the other, it displays a central double 
structure (as in Fig. 3d, left, for I odd). At lower T, 
however, only the double peak structure (features "C") 
remains (see Fig. 2c). The contribution from segments 
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with I even is seen to decrease rapidly with T, and one 
is led to conclude that a periodic structure with I odd 
becomes energetically preferable. 

A complete description of the NMR line at low T can 
now be proposed. A chain of N spins is considered, which 
contains a number 77 of bond-defects forming a periodic 
LAM structure with I odd. The dynamical behavior, 
i.e., the function K^v) in Eq. 2, is now to be considered 
as describing the dynamics (or the thermal fluctuations) 
of that whole periodic LAM structure To account 

for the additional static local fields induced by the chain 
edges, expressions for /i„ and /iAr_„ [given by Eq. (1)] 
are added to the applied field H in Eq. 2. The total 
NMR line, / = S„ /„, is then calculated by varying the 
parameters a, G, I and Dq. The parameters a and G 
are used to adjust the width of feature "B" . The shape 
of feature "C" is then determined by varying I and Dq, 
keeping in mind that the narrowing effect, i.e., the pa- 
rameter G, plays here a crucial role on the occurrence or 
not of such a small splitting at the centre of the NMR 
line. Finally the number ij of bond-defects (with the con- 
dition rjl < N) is used to adjust the intensity of features 
"A" relatively to that of the main line. With that pro- 
cedure, one obtains the open dots in Figs. 1 and 2. The 
agreement with the experiments is in general excellent. 
Particular confidence in this approach is provided by the 
high sensitivity of the final agreement to only a few pa- 
rameters, essentially I and G (compare the different lines 
in Fig. 3d, left). At low temperature, when features "C" 
are clearly visible, one always obtains the following de- 
termination: rjl ~ N {~ 1800 ± 200). This means that 
the periodic LAM structure develops practically over all 
the chain length. In that case, an evaluation of the den- 
sity of the bond-defects is simply given by p ~ For 
H / / SL, one obtains a rather accurate determination of 
the parameters a and G. They do not depend appre- 
ciably on T and H: a ^ 0.60 ± 0.04, G = (2.8 ± 0.2) 
10^"* GHz. For H / / c, the determination is much more 
uncertain (a ~ 0.46 and G ~ 4.2 10^'* GHz) as a small 
contribution with I even is ignored in our description . 
The parameter Dq remains also constant {Dq ~ 1.5 10^^ 
GHz for i7 // a, and Dq ~ 5.6 10"^ GHz for // c), 
except at the lowest temperature (~ 20 K) where larger 
values (by a factor ~ 2) provide a better agreement. At 
low T, however, the description represented in Fig. 3c 
becomes insufficient when the overlap of the LAM be- 
comes very large. The main variation is finally observed 
on the density of bond-defects p ~ /~^, which decreases 
almost linearly with T. This result is well established for 
H / / a (Fig. 2d) . The values obtained ioi H // c are 
also shown in that figure, though they are more uncer- 
tain for the reason explained above ||l^ . A first question 
raised by the present description concerns the source of 
the bond-defects (J' < J). Are they related to the inter- 
stitial excess oxygen which characterizes that compound? 
There are two reasons to rule out that possibility. First, 



at the lowest temperature, T = 20 K, the density of bond- 
defects (p ~ 2.8 10~^) remains more than ten times larger 
than the concentration of the residual spin- 1/2 impuri- 
ties (~ 1.3 10~^) after annealing. Second, the observed T 
dependence of p tells us that the number of bond-defects 
is not a constant but it decreases with T. Remarkably, p 
is observed to decrease linearly with T, as does the in- 
verse of the spin correlation length ~ 2T/J. That 
last remark strongly supports our assumption that the 
bond-defects are an intrinsic property of the "real" spin 
chains in Sr2Cu03. 

In the spin-Peierls phenomenology - and a priori above 
Tjv, Sr2Cu03 is a SPS in its U phase - the descrip- 
tion of the lattice is usually considered in its linear ap- 
proximation. The lattice modes are described by sim- 
ple phonon branches (characteristic frequency Vp) and 
the spin-lattice coupling is considered in its adiabatic 
limit [vp < J) [||. A more realistic description, however, 
should take into account the non-linear properties of the 
lattice. They may give rise to additional "localized" (and 
rather mobile) excitations (lattice solitons, for instance) . 
In the case of uniform chains, one kind of such non-linear 
excitations are of particular interest: the so-called "local- 
ized modes" as defined in [|l5|. They correspond to local 
oscillations of atoms, at high frequency and with large 
amplitude (breather like) . Their presence in the lattice of 
an antiferromagnetic chain would result in bond-defects 
as those discussed above |l^. Alternatively, if the SP 
phenomenology is considered in its non-adiabatic limit 
{vp > J), non-linearities might be also expected from 
the spin-lattice coupling itself. In that discussed 
recently, the SP transition, instead of being driven by a 
phonon softening, would be associated with a "critical" 
central peak ||l^, which usually signals the presence of 
non-linear excitations. Well above Tsp, could this effec- 
tive non-linearity explain the presence of localized exci- 
tations? The present NMR results may provide the first 
evidence of magneto-elastic localized modes in UQHC. 
A recent numerical investigation of such bond-defects 
in UQHC jl^ has already established the following re- 
sults: i) the expected magnetic interactions between two 
(static) bond-defects are attractive, ii) a LAM structure 
with I odd is energetically more favorable, iii) the LAM 
amplitude of an l-even structure decreases rapidly at low 
temperature. All these results agree with our analysis 
of the NMR line. Together with the possible repulsive 
interactions due to the lattice, they reinforce our propo- 
sition that a periodic LAM structure takes place at low 
T in UQHC. At that point, it is worth reminding that a 
double-peak structure, similar to our features "C", has 
been observed in the SPS CuGeOs well above Tsp, 
i.e., in its U phase. Our conclusion is that the presence of 
a LAM structure is strongly supported by the NMR re- 
sults. We propose a model - bond-defects and a periodic 
LAM structure - which explains both features "B" and 
"C" in an unique consistent model. That model deserves 



3 



certainly to be comforted experimentally and theoreti- 
cally. For instance, the quasi-static periodic AF structure 
of the LAM could be observed by neutron diffraction. 
Any description of such a "real" UQHC should consider 
explicitly the subtle relations between the spins and the 
lattice. 
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FIG. 1. NMR signals (solid lines) obtained at different fre- 
quencies u'^,ioi: H / / a compared to calculated values (open 
dots). "A", "B" and "C" denote the "features" to be ex- 
plained. In b), the dotted line describing feature "A" corre- 
sponds to = 0, while the calculation is obtained with (p = ^ 
(see dotted and solid lines in Fig. 3a). 



FIG. 2. a) NMR line for H // c. The inset shows the details 
of the central parts of the three hyperfine components. The 
open dots are calculated for a periodic LAM with / odd only, 
b) Example of a spectrum resulting from LAM with both / 
even (dotted) and odd (dashed line), with relative intensity 
0.13. c) Comparison between calculation for / odd only (open 
dots) and experiment (solid line) obtained at T = 20 K. d) T 
dependence of the defect density p for a model of a periodic 
LAM structure with / odd. The dotted line is a linear fit with 
the data for H // a. The data for H // c are also shown (see 
text). 



FIG. 3. a) Contour of the local field (in absolute value) as 
a function of the site number n apart from a bond defect: 
Y — h'y^/ {BAfsH) and ^, a and (j> are defined in text, b) De- 
scription of feature "B" by independent mobile bond-defects 
in the diffusive and ballistic models (lines with open dots) 
compared to experimental data (solid line), c) Contour of 
the local field between bond-defects with I odd and even as 
used in our calculations, d) NMR lineshapes resulting from 
periodic LAM with I odd and even: the solid lines are calcu- 
lated for the experimental conditions of Fig. lb, with G — 2.8 
10"" GHz for I = 241 (left) and I = 240 (right); the dotted 
(G = 2.2 10"* GHz, same I = 241) and dashed (/ = 261, 
same G = 2.810~* GHz) lines illustrate the high sensitivity 
of these parameters on the final lineshape. 
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